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SUMMARY 

A model for the analysis of vortex generators In a fully viscous subsonic 
Internal flow Is evaluated. A vortlclty source term Is used In a modified 
form of the "Parabolized" Navler-Stokes equations to model the shed vortex. 
Computed results are compared with Idealized flow vortex paths, and with 
experimental data for vortex generators embedded In a thick turbulent boundary 
layer. The analysis Is also compared with experimental data for a separated 
diffusing S-duct and for a diffusing S-duct with vortex generators. Quantlta 
tlve comparisons are shown for the latter three cases. Emphasis Is placed on 
verifying the ability of the model to predict global distortions In the flow 
field. 


INTRODUCTION 

Three-dimensional subsonic diffusers are very common In modern aircraft 
propulsion systems. Strong secondary flows are generated by the presence of 
offset bends and changes In cross-sectional shape of the duct. In addition to 
this phenomena, since the flow Is diffusing, the boundary layers may grow to a 
thickness comparable to the duct radius. These effects weigh heavily on the 
overall performance of the diffuser. 

In efforts to save weight and thereby fuel. It Is common to design these 
ducts to be as short as possible. In these cases the designer must be con- 
cerned with the flow separating due to a strong adverse pressure gradient. 

This separation may have devastating effects on engine performance, such as 
the engine unstarting. To alleviate this problem, vortex generators can be 
used as a flow control device. The vortex generator will cause a mixing of 
the higti momentum core flow with the low momentum flow In the boundary layer, 
resulting in a net increase of energy near the surface. This can delay or 
eliminate the separation depending on the flow conditions. There Is a penalty 
that Is paid for the use of the vortex generators and that Is the drag of the 
generator. But, since the vortex generators tend to have small planform 
areas, this adverse effect Is small. An excellent background on forced mixing 
can be obtained from references 1 to 3. 

To solve for the flow fields In these ducts, a complete three-dimensional 
Navler-Stokes analysis could be used, but large amounts of computer time and 
storage are required for a detailed solution. However, by making certain 
approximations to the Navler-Stokes equations, such as neglecting streamwise 
diffusion, a set of equations can be obtained for fully viscous internal flow 
that can be solved by forward marching In space. A computer analysis, desig- 
nated PEPSIG, has been developed using these equations to compute three- 
dimensional subsonic viscous flow through curved ducts with superel 1 Iptlc 



are both unity.) The known three-dimensional pressure field Is represented 
by pi = pi(X],X 2 ,X 3 ) and the unknown one-dimensional viscous blockage correc- 
tion Is p v = Pv( x 3 ) - 1° this approach, pi was computed using a three- 

dimensional potential flow analysis. This brings the elliptic effects of the 
geometry Into the solution. 

The secondary velocities v s and w s are split Into Irrotatlonal and 
rotational components. A scalar potential, <p, corresponding to the Irrota- 
tlonal component of the secondary velocity, Is defined as 
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The transverse momentum equations are redefined In terms of a streamwlse 
vortlcl ty-stream function formulation given by 
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Here Is the stream function, or vector potential, corresponding to the 
rotational component of the secondary velocity, and S2 Is the vortlclty In 
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The effect of the drag of the vortex generator Is also Included In the 
model. The drag on the wing section Is a combination of profile drag, which 
Is due to viscous and pressure effects, and Induced drag, which Is due to the 
shed vortex. In this model the profile drag of the vortex generator Is 
neglected In comparison to the Induced drag because In the cases studied here 
the generators were small. The Induced drag Is then proportional to the vor- 
tex strength and the crossflow velocity at a point. This term Is Included In 
the governing equations as a negative source term In the primary momentum 
equation. 


IDEAL FLOW RESULTS 

The first test case used to verify the vortex generator model was a com- 
parison of the motion of the center of a vortex due to the Induced velocity of 
another vortex. The results of the present analysis were compared with the 
results of an Ideal flow calculation with a wall simulated using Image vor- 
tlcles as shown In figure 3. In this method the velocity of a vortex, say 
vortex A In figure 3, can be calculated at any time by 



After the velocities at points A and B are calculated new positions can be 
determined by time Integration. In the viscous code this Ideal flow problem 
Is approximated by laminar Incompressible flow In a circular pipe. The 
Reynolds number based on pipe diameter was 2000. The boundary layers were 
thin, only 5 percent of the duct radius, with the vortex generator tip located 
well outside of the boundary layer. The total drag of the vortex generator 
was neglected. Four different cases are considered: (1) counter rotating 

vortlcles of the same strength that rotate toward each other with L/H = 2.5, 
(2) counter rotating vortlcles of the same strength that rotate away from each 
other with L/H = 2.5, (3) co-rotatlng vortlcles* of the same strength with 
L/H = 2.5, and (4) counter rotating vortlcles, one vortex having double the 
strength of the other, that rotate away from each other with L/H = 4.0. 

Figure 4(a) shows a comparison between the Ideal flow calculation of vor- 
tex paths and the calculated values of streamwlse vortlclty from PEPSIG for 
case 1. The vortlcles are rotating away from each other. The graph of the 
Idealized vortex motion shows the centers being attracted to each other for a 
short time and then they proceed to march away from the wall. The streamwlse 
vortlclty contour plots show the same type of motion of the vortlcles. Since 
the two vortlcles are moving toward the center of the duct the curved walls do 
not affect their motion. Although the viscous calculation was set up to 
attempt to model the Ideal flow the results of figure 4(a) show that the 
vortex strength Is dissipating a significant amount as It moves downstream. 
Figure 4(b) shows the same comparison for two vortlcles that are rotating 
toward each other. The graph of the Idealized vortex motion shows the centers 
being forced towards the wall, then the wall starts to force the two vortlcles 
away from each other. The four streamwlse vortlclty contour plots show the 
progression of the two vortlcles as they march down the duct. The vortlcles 
exhibit the same basic behavior as the Ideal flow results except for the con- 
fining effect of the curved walls. Again In these results the strength of the 
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agree well with the experiment. From this comparison It Is evident that the 
overall effects of the vortex are predicted but If detailed Information about 
skin friction and heat transfer are needed the effect of the wall on the vor- 
tex must be better modeled. 


DIFFUSING S-BEND WITHOUT VORIEX GENERATORS 

Figure 7 shows the 30 degree-30 degree S-bend diffuser tested experimen- 
tally In reference 13. The flow In this duct was turbulent with a Mach number 
of 0.6 and a Reynolds number based on the duct diameter of 1 760 400. The 
Initial conditions were measured at 1.65 duct diameters upstream of the first 
bend to remove the Influence of the bend on the static pressure. The Initial 
boundary layer thickness was 0.1 times the Initial duct radius and the area 
ratio was 1.51. In both the experiment and the analysis the flow In the duct 
separated. This Is due to the adverse pressure gradient of the area change 
combined with the effect of the pressure driven secondary flow due to the 
centerline curvature. The analysis marches through the separated region by 
using the "flare" approximation. A good explanation of the "flare" approxlma 
tlon Is presented In reference 14. 

Figure 8(a) shows the computed contours of a total pressure coefficient, 
Cpt = (Pt - Pa>)/q<x>* at various locations In the duct. Figure 8(b) shows the 
experimental contours of Cpj.. The maximum and minimum values of Cn^. at each 
streamwlse station are shown on figures 8(a) and (b). Comparing tne computed 
and experimental results Indicates that the analysis Is able to adequately 
predict the total pressure distortions for the duct. The reader must be aware 
that although the computed results In the separated region will not be correct 
because of the flare approximation, the global effect of the separated region 
Is well modeled. The comparison also shows that at 0 ■= 30 degrees, the 
fourth contour plot, the experimental results Indicate a larger separated 
region. In the experiment the streamwise separated region was found to be 
between e = 22 degrees and e = 44 degrees, while In the computation the 
separated region was between o = 30 degrees and 0 = 54 degrees. Figures 9 
and 10 show comparisons between the computed and experimental secondary veloc- 
ity profiles at the Inflection plane and at the bend exit. The comparison 
between the two results Is very good. At the Inflection point the vortex due 
to the curvature of the centerline Is evident In both plots. Also at the 
Inside of the first bend the separated region can be seen in the experimental 
results, by the region with no data. In the computed results the onset of 
separation Is also evident where there Is minimal secondary flow. At the bend 
exit both results from figure 10 show that the separated region Is gone by the 
large amount of flow being swept toward the outside of the second bend. These 
results differ from those of a nondiffusing circular cross section S-bend, 
where the vortex due to the centerline curvature Is strengthened In the second 
bend . 


DIFFUSING S BEND WITH VORIEX GENERATORS 

After the separation was found experimentally In the S-bend diffuser, 
three pairs of vortex generators were placed in the duct to remove the separa 
tion. Figure 11 shows the specifications for the vortex generators and 
figure 12 shows a schematic of the diffuser with the axial location of the 
vortex generators Indicated. The generators were placed well upstream of the 
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600 to 700 grid points per second. This Is approximately 150 to 200 times 
faster than a full Navler-Stokes analysis for comparable accuracy. 
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Figure 1. - A typical vortex generator. 
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(a) Co-rotating. 





(b) Counter rotating. 

Figure 2. - Typical vortex generator configurations. 
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(b) Case 2: Counter rotating vorticies of equal strengths rotating 
towards each other. L/H - 2. 5. 


Figure 4. - Continued. 
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(d) Case 4: Counter rotating vorticies with the strength of vortex 1 
twice that of vortex 2. Rotation is away from each other. L/H - 4. 0. 

Figure 4. - Concluded. 
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Figure 5. - Schematic of vortex generator locations in a straight duct. 
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Figure 10. - Comparison of secondary flow at the exit plane for S-duct without v.g.'s. 



(a) Dimensions of vortex generator. 
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(b) Dimensions of v.g. configuration. 
All dimensions in cm. 
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Figure 11. - Vortex generator configuration 
for S-duct diffuser. 




(b) Experimental total pressure coefficient contours for S-duct with v.g.'s. 
Figure 13. - Concluded. 



EXPERIMENTAL RESULT COMPUTED RESULT 


Figure 14. - Comparison of secondary flow at the inflection plane for S-duct with v.g.'s. 
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